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Abstract 



(N 

K*" ! We present a calculation of the 0{ma'^ In a) corrections to positronium 



energy levels. The result is used to estimate the current uncertainty in theo- 
fSJ ' retical predictions of the positronium spectrum. 

g ■ PACS: 36.10.Dr, 06.20.Jr, 12.20.Ds, 31.30.Jv 

^ I 1. Recently, a calculation of the positronium spectrum with 0{ma^) accuracy was 

,-J^ I completed [^0. The remaining uncertainty in theoretical predictions was estimated in 

I ■ using the value of the (9(ma^ln a) leading logarithmic corrections to positronium energy 

^' levels: 

■ ■ ' rr. ^499 ^ Nma^ln^a, ,^, 

XI 

^ I The purpose of this Letter is to present a detailed derivation of this result. 

In general, an appearance of logarithms of the fine structure constant in QED bound 
state problems is related to the fact that several momentum scales are involved in bound 
state calculations. Contributions that are logarithmic in a, usually appear as integrals of 
the form: 

^'^ F{k), (2) 



(2vr) 



3 



with F{k) r^ k ^ for the values of k such that ma <^ k <^ m ot ma^ <^ k <^ ma [§. 
Given the inequality k <^ m, it is possible to determine the leading logarithmic corrections 
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in the framework of the nonrelativistic Quantum Mechanics by expanding all perturbations 
in series of k/m and by using the time-independent perturbation theory. 

In order to find 0{ma^ \v? a) corrections, we first calculate all effective operators which 
deliver 0{ma'^\jia) contributions for n < 7. Such contributions first appear for n = 5 
and n = 6. Then, the (9(ma^ln a) corrections to energy levels are found by calculating 
relativistic corrections to these lower order operators. Some new operators arising at 0{ma'^) 
order should be also taken into account. For completeness we consider here the bound state 
of two different particles with masses m and M. 

2. Let us first consider the self energy operator of one particle in the Coulomb field of 
the other. The scattering amplitude^ reads: 



(27r)3 2A;^''^''^ ' k + H - E" 

The leading order logarithmic contribution is obtained by expanding the denominator {k + 
H — E)~^ in [H — E)/k and by using the leading non- relativistic approximation for the 
currents jj. Only spin-independent part of the currents should be considered since the spin- 
dependent part contains additional powers of k which destroy the logarithmic integration. 
One obtains: 

It is easy to see that the integration over k in the above equation is already logarithmic. For 
this reason we can neglect the non-commutativity of the Hamiltonian H with the exponent 
exp(iA;re). Fusing two exponents together, we integrate over k cutting the integral by the 
reduced particle mass /x = mM/{M + m) from above and by the particle energy ~ P^/a* 
from below. Including the self energy of the second particle, we obtain an effective operator: 



^(l-2J^]p{H-E)p\n^^. 
Snfi^ \ niM J ^^ 



y^^iP^r) -. --— l-2£-\v{H-E)v\^^- (5) 



Here and below we consider ln(p^/yU^) as commuting with all other quantities. 

By averaging this operator over nS states, one obtains the well-known non-recoil loga- 
rithmic contribution to the Lamb shift of the nS levels: 

^^ 3/i2 V mMj ^ ' 

In order to obtain 0{ma'^ \v? a) corrections induced by the operator Vg^^ it is necessary 
to calculate relativistic corrections to the expectation value of this operator. There are 
several sources of such corrections; below we analyze them. 

The simplest one is the relativistic correction to the currents. It is obtained by the 
substitution 



^The notations of [|[ are used throughout this paper. 
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performed for one of the currents in Eq.(Q). Again, only the orbital part of the currents 
contributes. The correction to the energy levels reads: 
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At the final stage of the calculation the self energy of the second particle was included. 
The next-to-next-to-leading order effect of the retardation is described by the operator 
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Though there is the fifth power of k in the denominator, the spin-dependent part of the 
currents still does not contribute to the double logarithmic correction. It can be easily seen, 
that the commutator [if, a x fcexp(ifcre)] is 0{k'^) as A; ^ 0, so that the resulting integral 
over k for the spin-dependent part of the currents is non-logarithmic. We therefore calculate 
the commutators in Eq. (^ and keep there only such terms that are quadratic in k. We 
obtain: 
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Next, we consider relativistic corrections to the Coulomb interaction and to the dispersion 
law of the particles in the intermediate state: 
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Taking the average value of this operator, one sees that the double logarithmic contribution 
is absent. 

Accounting for additional magnetic exchange between the particles requires some care. 
There exist eight irreducible diagrams. Only four of them, shown in Fig.l, deliver double 
logarithmic contributions. One finds, that the double logarithmic contributions from the 
diagrams Fig. la and Fig. lb compensate each other, while that of Fig.lc and Fig. Id sum up 
to the following energy shift: 
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Integration over k is performed in the limits (/ir^)"^ < k < fi. 

Finally, there exists 0{a'^) correction to the wave function of the bound state, i.e. the 
correction due to the iteration of the Breit Hamiltonian and the LO operator V^{p, r) (cf. 
Eq.® ). The calculation of this correction is described in the final part of this Letter, where 
all lowest order logarithmic operators are considered simultaneously. 

3. We consider now an exchange of a single magnetic photon between the two particles. 
The corresponding scattering amplitude is similar to Eq. (j^): 

- ^- = «/ (^^'^^(-^ ' -P' - g)^"^' % + ^^i, ^i(P - ^Py"""' + h-- (13) 

We used this amplitude in |0| when discussed the retardation effects. Considering the next- 
to-leading order retardation, one finds the operator 

where again the non-commutativity of e*^^ and H has been neglected. Because of the 
presence of the exp{iqr), the logarithmic integral over q in Eq.(|l^ is cut at 1/r ~ p from 
above 0. We arrive at the following effective operator: 



^'°^" -^ = -7^ -p{H - E)p\nP-^. 



Vlf'^p.r) = p(H-E)p\n^. (15) 



The sources of double logarithmic corrections to the single magnetic exchange are the 
same as in the case of the self energy operator. The only essential difference in two calcula- 
tions consists in a change of the upper cut-off. For this reason we skip a detailed discussion 
and present the results of the calculation. We obtain: 

^ 3mM \ mMj ' ^ ^ 

4. The next source of the double logarithmic corrections is the double magnetic exchange 
with two seagull vertices. The corresponding potential reads 



a^ f d^k in An ,„r, iee'^'^ 



^^ mM J {2nf2k2k' k + k' + H ~ E ' ^ ^ 

Here k' = q — k, e^ej = 6ij — k^kj/k"^ and similarly for e'. It is sufficient to consider the 
leading nonrelativistic approximation for the seagull vertex. From Eg . ([I9|) one sees, that 
the logarithmic contribution comes from the region of momenta where \q\ <S |fc| -C fi. 
Neglecting H — E in comparison with k and integrating over k from q to n one obtains: 

^L0^^1^i_ (20) 

We consider now relativistic corrections to the operator in Eq.(^). In this case there 
are other sources of corrections, as compared to two cases considered above. The following 
corrections should be considered - relativistic corrections to the seagull vertex, an appear- 
ance of the magnetic-magnetic-Coulomb vertex, an expansion of the "heavy" intermediate 
energy denominators in powers of |p|/wi, k/m, k'/m and the usuaf] retardation effects. We 
find that none of these effects produces the double logarithmic correction. 

For example, relativistic corrections to the seagull vertex are polynomial in the momenta 
p, k and q. On the other hand, the structure of the denominator in Eq.(|19[) does not change. 
In this case the logarithmic contribution comes only from the region of large |fc|, so that the 
resulting operator is of the form of Eq.(|20|) times a polynomial in external momenta. Such 
operators are too singular to produce a double logarithmic contribution. 

Hence we conclude that the C(ma^ln^a) corrections to the expectation value of the 
operator Eq.(^) appear only due to relativistic corrections to the wave functions. As we 
have mentioned already, these corrections are considered at the end of this Letter. 

5. We now turn to logarithmic operators which arise at the 0{ma^) order. The 0{a) 
correction to such operators can be caused only by the effect of retardation. There are two 
effects of this sort. Let us denote by k the momentum of the photon in the intermediate state. 
Then, for the "light" e~^e~'j intermediate state the energy denominators are expanded in 
{H — E)/k, whereas for the "heavy" e~^e~e~^e~'j intermediate state the expansion parameter 
is k/m. In contrast to these, all other relativistic corrections are of relative a^ order. 
Retardation corrections clearly require that the corresponding effective operators are induced 
by an exchange of at least one magnetic photon. Inspecting such operators in , we observe 
that only two of them, the "one-loop" operator and the operator induced by a double 
magnetic exchange with one seagull vertex, can give rise to (9(ma^ln^ a) corrections on one 
hand, and were not considered yet on the other. 

Retardation correction to the one-loop operator cannot produce the double logarith- 
mic contribution. Indeed, both effective vertices entering the one-loop operator, are p- 
independent, so that their commutators with H produce only positive powers of momenta. 
One can easily check that resulting operators are too singular to produce two logarithms. 

For the same reason considering the double magnetic exchange with one seagull vertex 
we keep orbital currents only. The sum of three diagrams (see Fig. 2) then gives: 



^By "usual" we mean here such retardation effects, that do not resolve the point-hke seagull 
vertex. 
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where 



^^^^^.fc'n^ ,^ = ,',,^qr.^ v, = ^e^^^^. (22) 

The contribution of three diagrams with the seagull vertex on the second particle line is 
obtained from Eq. (PTI) by the substitution m ^^ M. An account of the retardation to first 



order gives the operator: 

TA2ex a^ f d^k An A-n ( 1 1 rtr i , ^ ^ izj ^ 

+..^.3(^[//,..])+h.c.. (23) 

For k <^ q this operator reduces to 
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Double logarithmic correction to the energy arises if one integrates over k from g^/yU to g, 
and then over q from /ia to /i. The result is 

ir-E = -'-^^Wa. (25) 

^ 3mM ^ ^ 

The last "irreducible" correction arises due to the single-seagull diagrams with one of 
the magnetic quanta absorbed by the same charged particle (see Fig. 3). In contrast to the 
previous calculation, there is a double logarithmic contribution in the diagrams Fig.3a,b 
coming from the region of g -C /c -C /i. In the sum of Fig. 3a and Fig. 3b the contribution of 
this region cancels out. The result is : 

,-£ = _!£M5)!,„.„. (26) 

^ 3mM ^ ' 

Eq. (|26|) completes the analysis of the "irreducible" 0{ma^ \v? a) corrections. 
6. The last source of the (9(ma^ln a) corrections is related to the wave function modi- 
fication by relativistic effects: 

5^E = {VGU + UGV) . (27) 

Here G is the reduced nonrelativistic Green function, U is the Breit Hamiltonian projected 
on S'-states: 
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and the operator V is the sum of the lowest-order logarithmic operators Eqs.(^,|T^J2^ 
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The sum of the first two terms in U{p, r) can be represented as a linear combination of 
the operators H^, {H,C{r)}, and C(r)^ (see 0). The last two terms induce the following 
double logarithmic corrections: 
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In a similar manner we consider the contribution of the 6{r) operator from the Breit 
Hamiltonian. We obtain: 
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7. In order to check that we have accounted for all sources of the 0{m.a'^ In a) correc- 
tions, it is useful to consider the problem from a more formal point of view. As we mentioned 



above, logarithms of a arise because of the hierarchy of scales in QED bound states calcu- 
lations. In the framework of NRQED such logarithms may be detected by analyzing the 
on-shell scattering amplitude of two particles and identifying such contributions to this am- 
plitude that have zero index in the nonrelativistic region. Therefore, to get the 0{moi^ \y? a) 
corrections to energy levels, we have to consider the 0(a^) scattering amplitude and find all 
zero index graphs. 

To do so, we proceed in the following way: i) divide all graphs into classes according 
to the number of magnetic photons (e.g., class contains all graphs with four Coulomb 
exchanges, class 1 contains all graphs with three Coulomb and one magnetic photon, and 
so on); ii) in each class, we consider only such graphs that in the leading nonrelativistic 
approximation have negative or zero indices; iii) the negative indices are shifted to zero by 
accounting for retardation corrections as well as all possible relativistic corrections to the 
elements of a given graph. 

Consider first the class 0. In the leading non- relativistic approximation the only contri- 
bution comes from the four-Coulomb ladder graph with both electron and positron staying 
in positive-energy intermediate states. This graph has the index 3x3 — 3x2 — 4x2 = —5, 
where the first term comes from three-loop integration volume, the second one from three 
energy denominators, corresponding to e^e~ intermediate states and the last term arises 
from the four Coulomb propagators. To shift this index to zero, we have to consider rela- 
tivistic corrections with the fifth net power of momenta. However, all relativistic effects in 
pure Coulomb graphs bring in corrections which scale like 0[v'^\ i.e. they increase the index 
by 2. Therefore, we conclude that the class is free from potentially logarithmic graphs. 

Turning to other classes, we notice that one can significantly decrease the number of 
graphs by considering only such where all magnetic and seagull vertices belong to one and 
the same fermion line. It is easy to see that the equal-time transfer of a vertex from one 
fermion line to the other does not change the index of the graph. Hence, for the purpose 
of index counting we can consider the simplified problem of the O(a^) scattering of one 
particle at the external Coulomb field. After identifying potentially logarithmic graphs with 
magnetic and seagull vertices on one fermion line one must consider all possible equal-time 
transfers of those vertices between two fermion lines and all possible relativistic corrections. 

To make further discussion more transparent, we introduce the following notations. The 
Coulomb, magnetic, and seagull vertices are denoted by C, M, and S, respectively. For the 
class 1, in the leading nonrelativistic approximation we have the following graphs: MMCCC, 
CMMCC, CCMMC, and CCCMMQ. Their index, -3, is the sum of 9 from the integration 
volume, 2 from two Pauli currents in magnetic vertices, -6 from three Coulomb propagators, 
-6 from three energy denominators related to intermediate states without photon, -1 from 
the e"*"e~7 intermediate state, and -1 from the normalization factor of the propagator of 
the magnetic photon. To increase this index to zero it is necessary to account for the 
effects of retardation. Since each order of perturbation theory for the retardation increases 
the index by 1, we have to consider either the third order retardation effects or the first 



^The notations are such, that e.g. MMCCC represents a graph were first the magnetic photon is 
emitted, then it is absorbed and then three Coulomb exchanges occur. 



order retardation with additional (9(w^) relativistic corrections. Note that an account of 
the retardation includes not only an expansion of energy denominators {k + AE)~^, but 
also permutations of the Coulomb and magnetic vertices, e.g., MMCCC— >MCMCC. By 
inspection we find that all these effects were considered in our previous analysis of single- 
magnetic contributions. 

An analysis of the class 2 graphs is slightly more involved because of an appearance of the 
seagull vertex. The lowest index, -2, arises in the leading nonrelativistic approximation for 
the graphs with two seagull vertices, CCSS, CSSC, and SSCC. To compensate this index we 
have to account for the next-to- next-to-leading order retardation or some (9(f^) relativistic 
correction. The index -1 graphs include SMM, MSM, MMS, and MMMj^j with two 
additional Coulomb vertices, none of which is located either between S and M or Mj and 
Mj. The next-to-leading order retardation shifts the index of these graphs to zero. Finally, 
there are index graphs, MjMjMjMj and MjMjMjMj with two additional Coulomb photons, 
each located either to the left or to the right of all M-vertices. One can easily find the one- 
to-one correspondence between any of these graphs and some parts of the double-magnetic 
contributions considered in our previous analysis. 

In the class 3 there are only index graphs, SSMM and MMSS with one C outside pairs 
SS and MM. After equal-time transfer of one seagull and one magnetic vertex to the second 
fermion line, we get the graphic representation for the correction to the average value of the 
operator Eq. (pOD , induced by the magnetic part of the Breit Hamiltonian. 

Finally, an inspection of the class 4 graphs shows that all of them have positive indices. 

8. It remains to sum up all 0{ma'^\ii a) contributions. We obtain: 



5^avcr = :t— 4 - -r—v7 In «, (33) 




SEms = - ;7 1 - -^ In' a. (34) 

In the limiting case m/M — ^ 0, the old result for the hydrogen [^ is reproduced. Considering 
pure recoil corrections, one finds: 

It is interesting to note, that the 0{fi'^/{mMy) terms mutually canceled not only in Eq.(^), 
but also in individual contributions, i.e. those induced by the self energy, the single magnetic 
exchange and the double magnetic exchange with one seagull vertex. 

For the positronium, it is necessary to add the contribution caused by the virtual anni- 
hilation which is easily extracted from Eq. (p2|). One obtains: 



,3 + o-cr'aV(0) ,^2 
4 m'^ 



O^nnE = -3 — In a. (36) 



Using ^^(O)' = 5iom^a^/{87m^), we arrive at the final result for the 0{ma'^\n a) contri- 
bution to the positronium energy levels: 

^„ /499 „ ,\ ma^\Y?a^ ,„„. 

,B = - (_ + 7..' ^^^^^,0. (37) 



Its spin-dependent part is in agreement with the result of Ref . 0] . 

Numerically, the 0{ma'^\ii^a) contribution to the triplet energy levels equals to 

1 o 

6E (n^Si) = — J MHz. (38) 

The correction to the difference E{2^Si) — E{l^Si) amounts therefore to 1.16 MHz. For the 
singlet states, the 0{ma'^\n^ a) correction gives: 

5E (n^So) = -^ MHz. (39) 

We then find that the 0{ma'' In a) correction to the positronium ground state hyperfine 
splitting amounts to —0.9 MHz. 

We conclude, that the 0{ma'' In^ a) corrections turn out to be of the order of 1 MHz and 
hence somewhat enhanced, as compared to the naive estimate of the magnitude of the ma'^ 
effects. We note, that at 0{ma^) the remaining, non-logarithmic corrections, contributed 
approximately one half of the leading logarithmic ones. Extrapolating this situation to 
order 0{ma'^), we conclude that the current theoretical uncertainty in the predictions for 
the positronium energy levels |^,^ should be approximately 0.5 MHz. 
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Note added: We are indebted to K. Pachucki and S. Karshenboim for pointing out an 
error in Eq. (^) in the previous version of this Letter. After correcting this error, our result 
for the C(ma^log^ a) corrections to positronium energy levels agrees with that of Ref. M. 
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FIG. 1. Corrections to self-energy operator Eq.dsh due to additional magnetic exchange. 
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FIG. 2. Double magnetic exchange with one-seagull vertex. The dashed lines, representing 
Coulomb exchanges, show that the exact Green function for the system positronium plus photon(s) 
should be used in the calculation. 
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FIG. 3. Seagull corrections to single magnetic exchange. 
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